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The mass of the top quark (M t0 p) is a fundamen- 
tal parameter of the standard model (SM) and its 
large value makes the top quark contribution domi- 
nant in loop corrections to many observables, like the 
W boson mass Mw- Precise measurements of M\y 
and M top allow one to set indirect constraints on the 
mass of the, as yet unobserved, Higgs boson pj. 

In this Letter we present a measurement of M top 
using proton-antiproton collision events at a center- 
of-mass energy of 1.96 TeV. Top quarks are produced 
at the largest rate in pairs (tt), with each top quark 
decaying immediately into a W boson and a b quark 
nearly 100% of the timeQ. In this analysis events 
where both the W's decay to a quark-antiquark pair 
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(ti — > W + bW~b — ¥ qi^bq^qj)) are considered. This 
all-hadronic final state has the largest branching ra- 
tio among the possible decay channels (46%), but it 
is overwhelmed by the QCD multijet background pro- 
cesses, which surpass ti production by three orders of 
magnitude even after a dedicated trigger requirement. 
Nevertheless, it will be shown how this difficult back- 
ground can be successfully controlled and significantly 
suppressed with a properly optimized event selection. 
The fundamental analysis technique is the same ex- 
ploited to obtain the previous result from CDF, and 
is described in details inQ. However, improvements 
in the event selection and a larger dataset allow us to 
decrease the total uncertainty on M top by 21%. The 
additional dataset has been acquired at higher instan- 
taneous luminosity, which results in a higher number 
of background events in the data sample. Despite this 
fact, the introduction of significant improvements to 
the analysis results in the world best measurement of 
M top in the all-hadronic channel so far, also entering 
with the third largest weight in the M top world aver- 
age calculation (3,1^ ■ 

The data correspond to an integrated luminos- 
ity of 5.8fb _1 . They have been collected between 
March 2002 and February 2010 by the CDF detec- 
tor, a general-purpose apparatus designed to study 
pp collisions at the Tevatron and described in de- 
tail in Q . Events used in this measurement are se- 
lected by a multijet trigger 0, and retained only if 
they are well contained in the detector acceptance, 
have no well identified energetic electron or muon, 
and have a missing transverse energy $t satisfying 
t T l\JY. E T < 3GeV^, where £ E T is the sum of 
the transverse energy Er of all jets. Candidate events 
are also required to have from six to eight "tight" 
(E T > 15GeV and \r]\ < 2.0) jets. After this pres- 
election, a total of about 5.7 M events is observed in 
the data, with less than 9 thousand expected from 
ti events. To improve the signal-to-background ratio 
(S/B) a 6-tagging algorithm [8| is used to identify ("6- 
tag" or simply "tag") jets that most likely resulted 
from the fragmentation of a b quark. Only events 
with one to three tagged jets are then retained, im- 
proving the S/B by a factor of 6. In order to further 
increase the signal purity, a multivariate algorithm is 
implemented. An artificial neural network, based on a 
set of kinematic and jet shape variables [3|, is used to 
take advantage of the distinctive features of signal and 
background events. The neural network was trained 
using simulated ti events generated by Pythia [9| and 
propagated through the CDF detector simulation. At 
this level of selection the fraction of signal events is 
still negligible so that the data can be used to repre- 
sent the background. The value of the output node, 
-/V out , is used as a discriminant between signal and 
background, providing a gain in S/ B by an additional 
factor of about 30. 

The background for the ti multijet final state comes 
mainly from QCD production of heavy-quark pairs 



(bb and cc) and events with false tags from light- 
quark and gluon jets. Given the large theoretical 
uncertainties on the QCD multijet production cross 
section, the background prediction is obtained from 
the data themselves. The probability of tagging a 
jet in a background event (P + ) is evaluated using 
data with five tight jets and passing the preselection 
(S/B w 1/2000). This "tag rate" is parametrized in 
terms of a few relevant jet variables and is then used 
to estimate the probability that a candidate event be- 
longs to the background and contains a given number 
of tagged jets. As described in detail in Q this allows 
to predict the expected amount of background events 
in the selected samples as well as their distributions. 
For example, the average number of background 1-tag 
events is estimated by 



E 



Nu 



Ed 



tag 



n 
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where the outer sum runs over all events selected 
just before the 6-tagging requirement, and the in- 
ner one over the jets of the event. The factor Citag 
represents a correction to take into account correla- 
tions among jets within the same event and it is 
parametrized as a function of the same variables used 
for the tag rate. 

The analysis employs the template method to mea- 
sure M top with simultaneous calibration of the jet en- 
ergy scale ( JES) [1, [l(| , allowing a strong reduction of 
the associated systematic uncertainty. Distributions 
of variables sensitive to the "true" values of M top and 
JES, obtained by Monte Carlo (MC) events, are used 
as a reference ("template") in the measurement. A 
maximum likelihood fit is performed to define the val- 
ues that best reproduce the same distributions as ob- 
served in the data. An usual choice is to consider the 
distributions of the event-by-event reconstructed top 
quark mass, m™, and W boson mass, m^> c as the ref- 
erence templates. The JES is a multiplicative factor 
representing a correction applied to the raw energy of 
a reconstructed jet (E^™), so that its corrected energy 
E T = JES • i?^ aw , is a better estimate of the energy 
of the underlying parton [ll[ . Discrepancies between 
data and simulation result in an uncertainty on the 
JES value to be applied in MC events to reproduce 
the data, and, as a consequence, on the measurements 
of M top . Nevertheless, this value can be calibrated "in 
situ" , using m^ c as a template. This represents a well 
tested technique, first applied in [To| and now used to 
obtain the most precise top quark mass measurements 
at the Tevatron [3, [H| • 

The templates are built as follows Q . For each se- 
lected event, each of the six highest-i>r jets is assigned 
in turn to one of the six quarks of a ti all-hadronic fi- 
nal state. Then, for each combination the jets are 
arranged in two triplets (the top quarks) , each includ- 
ing a doublet (corresponding to the W boson) and a 
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b quark. To reduce the possible number of permu- 
tations, 6-tagged jets are assigned to b quarks only, 
resulting in 30, 6 or 18 permutations for events with 
one, two or three tagged jets, respectively [l2j . 

For each permutation m\ cc is obtained through a 
constrained fit based on the minimization of the fol- 
lowing x 2 -like function: 



Xt 



r 2 

1 w 



f , ( 



r 2 

mf/ b - m\ cc ) 



r? r 2 
2^ n* 



i=l 



where "^nff^ are the invariant masses of the two pairs 



(1 2) 

of jets assigned to light flavor quarks, frv--' b ' are the 
invariant masses of the triplets including one pair and 
one jet assigned to a b quark, Mw = 80.4GeV/c 2 
and Tw — 2.1 GeV are the measured mass and nat- 
ural width of the W boson [2|, and T t = 1.5 GeV is 
the assumed natural width of the top quark [l3| ■ The 
jet transverse momenta are constrained in the fit to 
the measured values, p™° as , within their known resolu- 
tions, Gi. The fit is performed with respect to m\ ec and 
the transverse momenta of the jets p^ i: and, among 
all the permutations, the one which gives the lowest 
value for the minimized x\ is selected. The variable 



m 



w 



is reconstructed by the same procedure consid- 
ered for mj ec , but with a \ 2 function, x\v> where also 
the W mass is left free to vary in the fit. The selected 



values of 



and 



enter the respective distribu- 



tions, built separately for events with exactly one or 
> 2 tags. 

Signal templates are built using MC events with 
M top values from 160 to 185GeV/c 2 , with steps of 
1.0GeV/c 2 , and, for each value, moving the JES by 
AJES-ojes from the default. Here ojes is the absolute 
uncertainty on the JES [ll| and AJES is a dimension- 
less number. Values of AJES between —2 and +2, in 
steps of 0.5, have been used, and in the following we 
refer to this parameter to denote variations of the JES. 
To construct the background templates we apply the 
fitting technique to the data events passing the neu- 
ral network selection cut, omitting the 6-tagging re- 
quirement ( "pretag" sample) Q . The weight of each 



value of 



and 



is given by the probability of 



the event to belong to the background and to con- 
tain tagged jets, evaluated by the tag rates of jets, as 
outlined above. 

Sets of simulated experiments ("pseudo- 
experiments" , PEs) have been performed to optimize 
the requirements on the values of V out , Xt an d Xw m 
order to minimize the statistical uncertainty on the 
Mtop measurement. As an improvement with respect 
toQ, two different sets of events, denoted by Sjes 
and SWtop > are used to build the m^> c and m\ cc tem- 



plates, respectively. The set Sjes is selected by using 
cuts on iV ou t and Xy/i while SM top is selected by a 
further requirement on Xt, so that SM top corresponds 
to a subset of Sjes- This new procedure contributes 
in reducingthe final total uncertainty on M top with 
respect to (3| by about 12%. TablesU and[ED report 
the flow of the event selection for 1-tag and > 2-tag 
events, respectively. As the final requirements are 
optimized separately for the two tagging categories, 
the 6-tag requirement is included in the flow just 
after the preselection. 



TABLE I: Selection flow for 1-tag events samples. For 
each requirement the number of events observed in the 
data, the expected number of ti signal events, the ab- 
solute efficiency on the signal (s) and the signal-to- 
background ratio (S/B) are shown. For the signal Mtop = 
172.5 GeV/c 2 and AJES = are used. The expecta- 
tions are normalized to the integrated luminosity of the 
data sample (5.8 fb" 1 ) using the theoretical cross section 
(7.46 pb), while the background is evaluated as the differ- 
ence between the data and the expected signal. 

Selection Requirement Data tt e(%) S/B 



Trigger + Presel. 
= 1 6-tag 
AU t > 0.97 

X 2 W < 2 (Sjes) 
X 2 t<HSM t J 



5683210 8854 20.6 1/641 

546579 3861 9.0 1/141 

5743 1028 2.4 1/4.6 

4368 881 2.1 1/4.0 

2256 604 1.4 1/2.7 



TABLE II: Selection flow for > 2-tag events samples. The 
same notations of Table|T] are used. 



Selection Requirement 


Data 


tt 


e(%) 


S/B 


Trigger + Presel. 


5683210 


8854 


20.6 


1/641 


> 2 6-tags 


47229 


1520 


3.5 


1/30 


A out > 0.94 


2379 


740 


1.7 


1/2.2 


Xw < 3 (Sjes) 


1196 


468 


1.1 


1/1.6 


X 2 t < 4 (S Mt J 


600 


316 


0.7 


1/0.9 



In order to measure M top with the simultane- 
ous calibration of the JES, a fit is performed in 
which an unbinned extended likelihood function is 
maximized to find the values of M top , AJES, the 
number of signal (n s ) and background (n&) events 
for each tagging category which best reproduce the 



observed distributions of 



and 



The 



likelihood depends on the probability density func- 
tions (p.d.f.'s) of m\ cc and m^- c expected for sig- 
nal (s) and background (&), P s {m\ cc \ M top , AJES), 
P s (m™ c \ M to p, AJES), A {mf% and P b (m^ c ). The 
notation points out that the shapes of the signal 
p.d.f.'s are functions of the fit parameters M top and 



G 



AJES. This dependence is obtained by fitting the 
whole set of templates, initially built as histograms. 
FigureQ] shows examples of signal and background 
templates for the > 2-tag sample, with the corre- 
sponding p.d.f.'s superimposed. 

tf m,'™ templates, > 2 tags events (AJES = 0.0) 
■ M lop = 160.0 
M lop = 172.5 
M lop = 185.0 
■P,(m"=|M AJES) 




250 300 
m,™ [GeV/c 2 ] 




£ — £l tag X £ 



>2 tags 



x £ A 



,TES C 



where £AjES co „ Bt r is a Gaussian term constraining the 
JES to the nominal value (i.e. AJES to 0) within its 
uncertainty : 



where, omitting the dependences on M top and Ajes, 

N S JES ___ roc ... rC c 



£ajes = ] I 

z=l 

N S ^o t 
ohs 

n 



n s P" lw (m w>i ) + n b P™ w (m Wti ) 



n b 



-M t0 p 



AsiisPT* (m t ,i) + AbiibP™* (m M ) 
A s n s + Abn b 



^ P(r a , n s ) ■ P(n, n b ) 



250 300 
[GeV/c 2 ] 



FIG. 1: Templates of m£ cc for events with > 2 tags and 
corresponding probability density functions superimposed. 
Top plot: the signal p.d.f, P s , for various values of M top 
and AJES = 0. Bottom plot: the background p.d.f., Pb- 

The presence of the different sets Sjes and SM top 
requires the generalizations of some of the terms of 
the likelihood with respect to [3J . The function can be 
divided into three parts: 



^B(< s , r„ A s ) ■ B (t b) r b , A b ) 



H< r b 



and 



£ 



N ' 



[ n b ™(b,BX P )J 
"(b.axp) 



In the first term the probability to observe the set 
m w, i, [i = 1) •••) N £f s ) OI m W va l" es reconstructed 
in the data is calculated by the signal and background 
expected distributions, P" lw and P™ w respectively, 
as a function of the free parameters of the fit M top , 
AJES, ri s , and n b . In the second the same is done 
for the distributions of the observed reconstructed top 

masses, mtj, (i = 1, —, ^obs'* " )> an d the m\ cc proba- 
bility density functions. The factors A s (M top) AJES) 
and Ab represent the acceptance of Su top with, re- 
spect to Sjes f° r signal and background, respectively 
(i.e., the fraction of events selected by the require- 
ments on Xt on ly)- F° r the signal this acceptance is 
parametrized as a function of the fit parameters M top 
and AJES. The third term, £ cv ts, gives the proba- 
bility to observe simultaneously the number of events 
selected in the data in the Sjes and the SM top sam- 
ples, given the assumed values for the average number 
of signal (n s ) and background (n b ) events to be ex- 
pected in Sjes and the acceptances ^4 s (M top , AJES) 
and Ab- It depends on the Poisson (P) and Binomial 
(B) probabilities 

G~ n ' Tl T 

P(r, n) = j — 



_ (JES-.TES conBtr ) 2 

£AJESeonstr = e 



[(jEScon 5 tr + AJES gJES )-JES constr ] 2 
2^ 

[AJES] 2 



Terms £itag and £>2tags are in turn defined as: 



£l,>2tags = £AJES x £M top X £cvts X C 



B (t, r, A) 



In the last term, C 



N° 



A t ■ (1 - A) r 



the parameter n b is con- 



N ' 



strained by a Gaussian to the a priori background 
estimate i.e. n^ b ,cxp) — 3652 ± 181 for 1-tag events 
and ri(b : oxp ) = 718 ± 14 for > 2-tag events. 

The possible presence of biases in the values re- 
turned by the likelihood fit has been investigated. 
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Pseudo-experiments are performed assuming specific 
values for M top and AJES and "pseudo-data" are 
therefore extracted from the corresponding signal and 
background templates. The results of these PEs have 
been compared to the input values, and calibration 
functions to be applied to the output from the fit have 
been defined in order to obtain, on average, a more re- 
liable estimate of the true values and uncertainties. 

Finally, the likelihood fit is applied to data. After 
the event selection described above, we are left with 
4368 and 1196 events with one and > 2 tags (147 have 
3 tags), respectively, in the Sjes sample. The corre- 
sponding expected backgrounds amount to 3652 ± 181 
and 718 ± 14 events, respectively. The tighter require- 
ments used for the SM tav samples select 2256 with one 
tag and 600 with > 2 tags (76 have 3 tags), with av- 
erage background estimates of 1712 ± 77 and 305 ± 22 
events, respectively. 

For these events the variables m^p and m\ ec have 
been reconstructed and used as the data inputs to the 
likelihood fit. Once the calibration procedure has been 
applied, the measurements of M top and AJES are 



M top 
AJES 



172.5 ± 1.4 (stat) ± 1.0 (JES) GeV/c 2 
-0.1 ± 0.3 (stat) ± 0.3 (M top ) • 



Figured] shows the measured values together with 
the negative log-likelihood contours whose projections 
correspond to one, two, and three a uncertainties on 
the values of M top and AJES as obtained from the 



likelihood fit. 



<o 1.5 



0.5 



-0.5 



-Ln(L/L ) Contours , > 1-tag events (5.8 fb" 1 ) 



>K Fitted Values 

-Ln(UL m=> ) = 4.5 

- -Ln(L/L m=> ) = 2.0 

-Ln(L/L mi> ) = 0.5 

i i i I i i i I 




-1.5 

166 168 170 172 174 176 178 180 

M,^ [GeV/c 2 ] 

FIG. 2: Negative log-likelihood contours for the likelihood 
fit performed for the Mt op and AJES measurement. The 
minimum is shown along with the contours whose projec- 
tions correspond to one, two, and three a uncertainties on 
the Mtop and AJES measurements. 

Figure [3] shows the m\ cc and m^p distributions for 
the data compared to the expected background and 
the signal for M t0 p and AJES corresponding to the 
measured values. The signal and background distribu- 
tions are normalized to the respective yields as fitted 
to the data, with the 1-tag and > 2-tag contributions 
summed together. 
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FIG. 3: Distributions of m r t cc (top plot) and m 
plot) as obtained in the selected data (black points) with 
> 1 tags, compared to the distributions from signal and 
background corresponding to the measured values of M top 
and AJES. The expected distributions are normalized to 
the best fit yields. 



Various sources of systematic uncertainties affect 
the M top and AJES measurements, as described in 
They are evaluated by performing PEs using tem- 
plates built by signal samples where effects due to 
systematic uncertainties have been included. The dif- 
ferences in the average values of Mt op and AJES with 
respect to the PEs performed with default templates 
are then taken into account. Possible residual biases 
existing after the calibration, and uncertainties on the 
parameters of the calibration functions are also taken 
into account. The largest contributions come from un- 
certainties on the modeling of the background, on the 
simulation of tt events, and on the individual correc- 
tions which JES depends on [ll[ . Table IIIII shows a 
summary of all the systematic uncertainties. 

In summary, we have presented a measurement of 
the top quark mass in the all-hadronic channel, using 
pp collision data corresponding to an integrated lumi- 
nosity of 5.8 fb -1 . An optimized event selection, based 
mainly on a neural network and a 6-tagging algorithm, 
allows us to select candidate event samples with S/B 
close to 1 in spite of the huge background still exist- 
ing at trigger level. The simultaneous calibration of 
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TABLE III: Sources of systematic uncertainty affecting 
the M top and AJES measurements. The total uncertainty 
is obtained by the quadrature sum of each contribution. 



Source 


6M top 5AJES 




(GeV/c 2 ) 



0.03 
0.01 
0.21 
0.04 
0.05 
0.01 

0.04 
0.04 
0.12 
0.05 
0.11 

01)4 

Total 1.1 0.29 



the jet energy scale, following a well established tech- 
nique, allows to reduce down to 1 GeV/c 2 the system- 
atic uncertainty due to this source. The value ob- 
tained for the JES is in agreement both with the de- 
fault value [ll[ and with the results obtained by other 



measurements of the top quark mass performed by the 
CDF Collaboration using the in situ calibration tech- 
nique [3,0]. The measured value of the top quark mass 
is M top = 172.5 ±1.4(stat) ±1.0 (JES) ± 1.1 (syst) 
GeV/c 2 , with a total uncertainty of 2.0GeV/c 2 . This 
result complements and is consistent with the most 
recent measurements obtained in other channels by 
the CDF and DO Collaborations, and also represents 
the most accurate all-hadronic measurement at the 
Tevatron so far. 
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